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Abstract 

Purdue and Favata calculate the tidal heating used certain classical pseu¬ 
dotensors. Booth and Creighton employed the quasi-local mass formalism of 
Brown and York to demonstrate the same subject. All of them give the result 
matched with the Newtonian theory. Here we present another Hamiltonian 
quasi-local boundary expressions and all give the same desired value. This 
indicates that the tidal heating is unique as Thorne predicted. Moreover, we 
discovered that the pseudo-tensor method and quasi-local method are funda¬ 
mentally different. 


1 Introduction 

In gravitation, one of the interested topics is that calculate the tidal heating: the 
interaction between a nearly isolated gravitating body and the external universe. 
Solar system provides a typical example for the illustration: Jupiter and its satellite 
Io [Tj|. Purdue [2] and Favata f3] examined the tidal heating for the classical pseudo¬ 
tensors. Booth and Creighton used the quasi-local mass formalism of Brown and 
York to demonstrate the same subject [4j. All of them give the same value as the 
Newtonian perspective. Here we present another Hamiltonian quasi-local boundary 
expressions to examine the tidal heating, we find that the result is unique as Purdue 
achieved, i.e., quasi-local expressions independent. 

Thorne claimed that all pseudo-tensors give the same tidal work as the Newtonian 
gravity |5]. Nester realized that pseudo-tensor method and quasi-local formalism are 
basically the same [6], however, Booth and Creighton prefer the quasi-local method 
such that all quantities can be manipulated in terms of real tensors on the quasi-local 
surface [3]. Thus one may imagine that there is no surprise using quasi-local expres¬ 
sions to give the same desired tidal work. Although using the pseudo-tensor method 
and quasi-local method give the same tidal heating rate, we find that the fundamental 
principle between these two are different. In particular, the Mqller pseudo-tensor give 
the standard tidal heating [3] but failed for inside matter requirement [7j . Meanwhile 
the Hamiltonian quasi-local method can rehabilitate this handicap. 

Generally speaking, different energy-momentum pseudo-tensor refers to different 
gauge condition for the gravitational energy localization, while different quasi-local 
boundary expression focus on different boundary condition. Confined to the tidal 
heating, we claim that the gauge condition and boundary condition are elementarily 
the same terminology. Different gauge condition corresponds different E int , where E int 
is the energy interaction between the isolated planet’s quadrupolar deformation and 
the external tidal held. We find that the tidal heating remains unchange for different 
quasi-local boundary expressions, thus the tidal heating is gauge invariant. However 
there is a change: the exchangeable energy rate E int , they are gauge dependent. 
Here we explain some terminology used in present paper. The expected tidal heating 
or tidal work rate W = —| IijE l \ where W refers to tidal work, the dot means 
differentiate w.r.t. time t, I t j is the mass quadrupole moment of the isolated planet 
and Eij is the tidal held of the external universe. Both A,- and Eij are time dependent, 
but symmetric and trace free. Here we emphasize that the tidal heating at which the 
external held does work on the isolated body and this is an energy dissipation process 
which means energy unexchangeable. Conversely, there is an exchangeable process 


E int ~ IijE u ). Purdue uses ^ nt = ^y£jE lJ to distinguish different options how 

to localize the gravitational energy by tuning the coefficient /3 | 2 ]. 

2 Technical background 

We used the same spacetime signature and notation as in [ 8 ]: let the geometrical 
units G = c = 1, where G and c are the Newtonian constant and speed of light. The 
Greek letters refer to the spacetime and Latin letters indicate the spatial. For the 
idea of energy-momentum pseudo-tensor tJ 1 , choose an appropriate super-potential 

u a w 

d„U a [H = 2 n^{T a » + £/), ( 1 ) 

where is the stress tensor, d v U a ^ can be described as the total energy-momentum 
complex and it is conserved since = 0. Recall the Einstein equation: 

G^ v = where G MI/ is the Einstein tensor and n = 87 rG/c 4 . Look closer at 

d consists two parts: 2 G a M is the mass energy inside matter and tj 1 ' is the 
gravitational energy-momentum in vacuum. The component t 0 J is the gravitational 
energy flux density. The criterion of the interior mass-energy is important. In partic¬ 
ular, the classical Mpller pseudo-tensor cannot satisfy this inside matter condition [7]. 
Thus one can conclude that Mqllcr pseudo-tensor is not appropriate to describe the 
energy-momentum in vacuum. 

The gravitational tidal heating rate can be computed as 

2kW = j y/-gt 0 j hj r 2 dQ, (2) 

where r = \J8 a bX a x b is the distance from the body in its local asymptotic rest frame 
and hj = Xj/r is the unit radial vector. In our calculation, the metric tensor can 
be decomposed as g a p = r] a p + h af 3 and its inverse g a > 3 = g a ^ — h a/3 . We have the 
following physical expressions [2]: 

2 M 3 2 • 10 • 4 . 

ho 0 = + -5 hbX a x b - E ab x a x b , h 0j = -p^hjX 1 - — EabX^Xj + —E^r 2 , (3) 

note that h t j = 8ijh 00 . 


3 Hamiltonian quasi-local boundary expressions 

Geometric theories are invariant under local diffeomorphism. Here we review the 
Hamiltonian quasi-local boundary expressions from a first order Lagrangian [9]: 

£ = dq Ap- A{q,p), (4) 

where q and p are canonical conjugate form fields, A is a potential. Let q be a /-form 
and e = (—1)A The corresponding Hamiltonian 3-form is defined as follows 

H(N) := £ N q Ap — (5) 

Taking the interior product of the Lagrangian density 

i N C = £ N q A p — ei N q A dp — edq A i N p — i N A — d(i N q A p), (6) 

where the Lie derivative £n = i^d + di^- Using (J 6 ]) , rewrite (l5j) 


'H(N) = + dB(N), 


( 7 ) 




where N ll dip = ei^qAdp^edqAi^p-Ai^A which is proportional to the held equations 
and vanishes ‘on shell’. Note that N M is the vector held. The Hamiltonian density 
din determines the evolution equations and initial value constraints. The natural 
boundary term B(N) = i^q A p. However, this boundary term is not unique since it 
can be removed by introducing a new Hamiltonian 

n'(N)=n(N) + d(-i N qAp) = N‘ i n ii . (8) 

Taking the variation of this new Hamiltonian 

8di'{N) = -ijv(F.E.) -Sq A £ N p + £ N q A Sp + dB(N), (9) 

where the held equation F.E. = Sq A A Sp. The boundary variation term is 

B(N) =—i N q A Sp + eSq Ai N p. (10) 

This B(N) cannot be removed because it comes from Sdi! directly. Boundary condi¬ 
tions can be obtained through the boundary term in the variation of the Hamiltonian 
vanishes. We add an appropriate boundary term to the Hamiltonian 

n'{N) -A H k (N) = N^din + dB k (N), ( 11 ) 

to modify the variational boundary term. In order to achieve nice components like 
iN^Sq A Ap) or i^(Aq ASp), there are four simple boundary expressions can be added. 
The variation of the four Hamiltonians including this four expressions are 

Sdi q (N) = K + diN(8q A Ap ), 5did{N) = K — d(zjvAg A Sp — eSq A *jyAp),(12) 
5di p (N) = K — di]y(Aq A dp), 5di c (N) = K + d{i^Sq A Ap — eAq A ijv<5p),(13) 

where K = — i N (F.E.) — 5q A £mP + £Nq A 8p. Thus we recovered ( 6 ] 

B q (N) = i N q A Ap — eAq A i N p, Bd{N) = i N q A Ap — eAq A i^p, (14) 
B P (N) — i N q A Ap — eAq AznP, B c (N) — i^q A Ap — eAq A InP, (15) 

where A q = q — q, Ap = p — p, both q and p are the background reference values. 
Alternatively, rewrite the above four equations in a compact form 

B kl ,k 2 ( N ) = Bp( N ) + kii N Aq A Ap + ek 2 Aq A i N Ap, (16) 

where k\ and k 2 can be 0 or 1. In detail Bop = B p , £>o,i = B d, B\p = B c and B\,\ = B q . 

Using the analogy of classical electrodynamics and apply to the relativistic gravity, 
the type of boundary condition should be either Dirichlet or Neumann, and even a 
mixture of these two. From the boundary condition point of view, we prefer B q and 
B p since they are the simplest, i.e., Dirichlet or Neumann. Meanwhile, the boundary 
conditions of B c or Bd could be a certain linear combination of Dirichlet and Neumann. 
For the case of B q , there are two ways to obtain the variation of the boundary term 
that satisfy iN^dq A Ap) = 0. First control q, then Sq = 0. The second is to freely 
vary q and then Sq becomes arbitrary, which implies Ap = 0. This is the ‘natural 
boundary condition’ as it forces p — p. Similarly for the variation boundary term 
iN(Aq A Sp) = 0. In addition, we have modified (ki,k 2 ) -A (ci,c 2 ), where Ci and c 2 
are arbitrary constants, such that the variation is still legitimate |Sj. 


3.1 Quasi-local M0ller and Freud super-potentials 

Here we apply this Hamiltonian formalism to the Einstein-Hilbert Lagrangian [6] 


£gr R a p A rj a ^ 


(17) 


where the curvature 2-form R a p = dY a p + r" 7 A T 7 ^, the connection 1-form Y a p = 
Y a p 1 dx y and the dual basis is rf-P"' = *(dx a Adx 13 ■ ■ ■). As before, the interior product 

i N C = i N R°p A r lo P + R a p A i N rjJ = £ N Y a p A -R(N). (18) 

After a straightforward manipulation, the Hamiltonian density from above becomes 

H(N) = N P U P - i N Y a pDgJ + dB(N ), (19) 


where R a = 2G p a g p which satisfies the dynamical evolution and initial value con¬ 
straints, i nY p u A Drfrj, 1 ' vanishes since the metric compatibility. Finally the boundary 
term B = A g p v = N a g ua Y p aa g pi/ . Note that it is legal to modify the boundary 

term replace a negative sign, i.e., B —> —B. Rewrite (TTilj) 


H(N) = N c 


2 G p aVp - -V^ 9 { 9 va Y\ a - g^Y\ a )dS l 


fils 


( 20 ) 


where dS pu = \e pv ^ K dx^ A dx K . Looking at (l20lh the term M U a ^ = 2yJ—gg a ^Y u ^ aa 
looks like Mpller super-potential [10] but not exactly. Instead we call this quasi-local 
Mollcr super-potential. The reason is that the Mpllcr pseudo-tensor cannot fulfill the 
interior mass-energy requirement 2 G p p as explained before, but the quasi-local Mpllcr 
super-potential uG a AA can. The variation of (U9|) 

5R = field equation terms + di^idY 01 p A g a ^). (21) 


This can be classified as Neumann type boundary condition, i.e., the connection 
T ~ dg is to be held fixed. Alternatively this terminology can be translated as the 
deDonder gauge: 0 = d K (^/—gg^ K ) = y/—gY^ K K . Using the analogy of the pseudo¬ 
tensor method, the quasi-local Mpller pseudo-tensor is MtaA = ^(mWq^). It is 
known that the tidal heating is Wm = — ^RjE 1 ^ where E - mt = ^^^( RjE u ) = 0 
which means the energy localization chosen (3 = —2 [3]. 

Keep the same track in [6j, swap the quasi-local Mpller super-potential to the 
other pattern B' = Y p u A i^r\ p = \N a gP a Y x p^l^g. Again flap the sign of B' and 
rewrite (HI 


n(N) = N a 


2 G p a r lp 




( 22 ) 


where pldj^ = —\f—gg^ a Y x p 1 d'l^ l ' a looks like the Freud super-potential [lTj but we 
call this the quasi-local Freud super-potential because it comes from the Hamiltonian 
formalism. The quasi-local Einstein pseudotensor is pt a p = d v (pU a ^). The known 
result for the tidal heating Wp = — , where F int = ^|(/yU J ) = 0 

which means the energy localization selected j3 = 1 [3]. Note that this boundary 
condition can be described as the Dirichlet type which means fixing \J—gg^ a - 

Based on the pseudo-tensor method, Thorne claimed that the tidal heating is 
unique [5] and Purdue verified that indeed it is gauge-invariant [2]. As far as the 
tidal heating is concerned, we find that the gauge condition and boundary condition 
are equivalent. Moreover we are going to verify that all the quasi-local boundary 
expressions obtain the standard tidal heating rate. 






3.2 Relativistic quasi-local boundary expressions 

Here we write the modified quasi-local expressions in holonomic frames [9] 


B(N) = B p (N) + Cl i N AT a p AAriJ-c 2 AT a p Ai N Ar 1 J 

= -l N a ( F Uj^+c ly /^h^T\JZ + C2^h^T\^Z)^ (23) 

where AT 0 ^ = r “^ Ci,c 2 are real and finite. For simplicity, consider the 

reference for flat spacetime F"^ = 0 in Cartesian coordinates. Extract the quasi-local 
superpotential in (f23ll 

= fW 0 h + Cl V^h x *r\JZ + C 2 y/=gh*r r w5lz (24) 


Bear in mind that the Hamiltonian has already fulfilled the inside matter value 2G a ■ 
This quasi-local Freud super-potential and the extra higher order terms hT only 
contribute the energy-momentum in vacuum. Carry on the calculation and we have 
the quasi-local pseudo-tensor 


t M 


E tZ + ci[(r A \ + T* x u )T° a v + h Xn r a a7rtU ]5Z 
+c 2 [(r^ + v a Z)B T \p + h^r T x ^\5 x z 

<$£(r*W - r\ v r'^) + v% a {v^\ - t"Z) + (r^ a + r^ a )r' vfi - 2r^ Q r% 
+cr [r A \r^ A/3 + tP\t v z - + hrv v omjV - h™ rv. 


+C 2 


SZ2T^ X F X ^ - V% X Y X \ - T% X T" X U ) + (2T^ a + T^ a )T\p 

-r^ a (3rv + V) + - h^R% av 


(25) 


Consequently the tidal heating is 


W 


3 + 2C2 Uh j E ij ) ~ 

10 dV 1 ’ 2 3 


(26) 


This shows that the term with c 2 contribute non-vanishing E- mt . In contrast, the 
term with c\ contributes nothing. Hence, the tidal heating rate is indeed boundary 
conditions independent. 


4 Conclusion 

Purdue and Favata calculate the tidal heating used classical pseudo-tensors. Booth 
and Creighton employed the quasi-local mass formalism of Brown and York to demon¬ 
strate the same subject. All of them give the result matched with the Newtonian 
gravity. Here we present another Hamiltonian quasi-local boundary expressions and 
all give the same desired result. This illustrates that the tidal heating is unambiguous 
since it is gauge invariant. In fact the ambiguity comes from E mt : energy interaction 
between the isolated body quadrupolar deformation and external tidal field. 

Meanwhile one may argue that this tidal heating uniqueness is natural because the 
principle of the pseudo-tensor method and quasi-local method are essentially equiva¬ 
lent. Although the main purpose of these two methods is to evaluate the tidal heating 
and obtain the same value, we discovered that they are different fundamentally. We 
prefer the Hamiltonian quasi-local formalism since it guarantees the dynamical evo¬ 
lution and initial value constraints, while the pseudo-tensor method cannot warranty 
the inside matter requirement. 






References 

[1] Peale S J, Cassen P and Reynolds R T 1979 Science 203 894 

[2] Purdue P 1999 Phys. Rev. D 60 104054 

[3] Favata M 2001 Phys. Rev. D 63 064013 

[4] Booth I S and Creighton J D E 2000 Phys. Rev. D 62 067503 

[5] Thorne K S 1998 Phys. Rev. D 58 124031 

[6] Chang C C, Nester J M and Chen C M 1999 Phys. Rev. Lett 83 1897 

[7] So L L, Nester J M and Chen H 2009 Class. Quantum Grav. 26 085004 

[8] Misner C W, Thorne K S and Wheeler J A 1973 Gravitation (San Francisco, CA: 
Freeman) 

[9] So L L 2007 Int. J. Mod. Phys. D 16 875 

[10] M0llcr C 1958 Ann. Phys. NY4 347 

[11] Freud Ph 1939 Ann. Math. 40 417 



